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ABSTRACT 



Context. LS I +61 303 is a member of the select group of gamma-ray binaries: galactic binary systems that contain a massive star 
and a compact object, show a changing milliarc second morphology and a similar broad spectral energy distribution (SED) that peaks 
at MeV-TeV energies and is modulated by the orbital motion. The nature of the compact object is unclear in LS I +61 303, LS 5039 
and HESS J0632+057, whereas PSR B 1259-63 harbours a 47.74 ms radio pulsar. 

Aims. A scenario in which a young pulsar wind interacts with the stellar wind has been proposed to explain the very high energy 
(VHE, E > 100 GeV) gamma-ray emission detected from LS I +61 303, although no pulses have been reported from this system at 
any wavelength. We aim to find evidence of the pulsar nature of the compact object. 
fG ' Methods. We performed phased array observations with the Giant Metrewave Radio Telescope (GMRT) at 1280 MHz centred at 

phase 0.54. Simultaneous data from the multi-bit phased array (PA) back-end with a sampling time of f samp = 128 >c/s and from the 
polarimeter (PMT) back-end with £ samp = 256 //s where taken. 

Results. No pulses have been found in the data set, with a minimum detectable mean flux density of ~ 0.38 mJy at 8-cr level for the 

■ pulsed emission from a putative pulsar with period P >2 ms and duty cycle D = 10% in the direction of LS I +61 303. 
Conclusions. The detection of posible radio pulsations will require deep and sensitive observations at frequencies ~ 0.5-5 GHz and 
orbital phases 0.6-0.7. However, it may be unfeasible to detect pulses if the putative pulsar is not beamed at the Earth or if there is a 
strong absorption within the binary system. 

Key words, binaries: close - gamma rays: stars - pulsars: general - radio continuum: stars - stars: individual: LS I +61 303 - X-rays: 

■ binaries 

irv 

f^i ' 1 . Introduction O rbital X-ray periodicity has also been fou nd dParedes et al.l 

, --J ! j : 1 Il997h lHarrison et all 120001 iTorres et aTll2010h . Similar results 

. Located at a distance of 2.0 ± 0.2 kpc (| Frail & Hjellming i have been obtained at hi gher energie s with INTEGRAL data 

^ ■ U221l>, LS I +61 303 contains a rapidly rotating BOVe star dHermsen & Ku ioer 2006; Zh ang et al.ll2010h . LS I +61 303 is 

with a stable equatorial shell, and a compact object of un- a i so spatially coincident with a high ene rgy (HE, E > 100 MeV) 

known nature with a mass between 1 and 5Mp, orbiting g amma -ray source detected by EGRET (KniffenetaDQWzl)- 

it every 26.5 days (Hutchings & Crampton 1981; Casares et al. _,._.„ ^ , T . _ , 

2005 4 Optical and IR orb ital modulatio n have been found , r Th * ^ rmi S P" ce u Telescope Large Area Telescope 

(Mendelson & Mazeh 1994; ParedesISD[l993). Using radial (Fermi/LAT) reported the first detection of the orbital mod- 

S velocities data, Casaresj*alJ (|20oSffound that the orbit is ec- ulatlon a J ™ w u lth u a P enod of 26 f>± M4m™ m 

centric (e - 0.72) and periastron takes place at pha se 0.23 + 0.02, anti-correlated with the X-ray emission (|Abdo et al. |l2009a|). 

taking T = JD 2, 443, 366.775. More recently lAragona et all At very high energy (VHE, E> 100 GeV) gamma rays, 

( 2009) stated that the eccentricity is slightly lower (e * 0.54) and LS I +61 303 was detected by the MAGIC Cherenkov telescope 

the phase of the periastron passage is 0.275 . The compact object dAlbert et al. 2006) and confirmed by the VERITAS stereoscopic 

would be a ne utron star for in clinations 25° < i < 60° and a black array (|Acciari et al.||2009a|). Further observations by the MAGIC 

hole if i < 25° dCasares et al.ll2005al) . However, ther e are discrep- collaboration led to the discovery of the orbit al variability ofTeV 

ancies b etw een the proposed orbi tal solutions bv ICasares etaD emission with a period of 26.8 + 0.2 days ([Albert et al.||2009|). 

(I2005ah and lAragona etalJ d2009h , and the inclination angle is Simultaneous observations of MAGIC and the XMM-Newton 

poorly constrained by radial velocities of the companion star and Swift X-ray satel lites revealed a correla tion between the X- 

alone. ra ys and VHE bands dAnderhub et al.l l2009). 

Quasi-periodic radio outbursts monitored during 23 years iMassi et aD (120041) reported the discovery of an extended jet- 

have prov ided an accur ate orbital period value of 26.4960 + like and precessing radio emitting structure at angular extensions 

0.0028 d (Gregory 2002). The maximum of the radio outbursts of 10-50 milliarcseconds. Owing to the presence of apparently 

varies between phase 0.45 and 0.95. relativistic radio emitting jets, LS I +61 303 was proposed to be a 
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microquasar. However, VLBA images obtained du ring a full or- 
bital cycle show a rotating elongated morphology ( Dhawan et al.l 
2006), which may be consistent with a model based on the in- 
teraction between the relativistic win d of a young non-accreting 
pulsar and th e wind of the donor star dMaraschi & Trev es 1981; 
iDubusI 12006b . On 2008 September 10, the Swift Burst Alert 
Telescope (BAT) detected a hard X-ray burst in the direction of 
LS I +61 303 which, assuming the association, would be the sig- 
nature of a magnetar-like activity that has been proposed to be 
linked to the pre sence of a young highly magnetized pulsar in the 
binary system (iDubus & Giebelsl [20081) . The Fermi/LAI spec- 
trum is compatible with a power law and an exponential cutoff 
at ~ 6 GeV, suggesting that there are two separate spectral com- 
ponents for the HE and VHE emission. Moreover, the spectral 
similarity with gamma-ray pulsars leads to the consideration of 
a magnetospheric origin for this HE component despite the fact 
that no pulses have been detected and although the orbital mod- 
ulation would b e unexpected in this scenario (lAbdo et al. 2009a; 
Petr i & Dubusl 1201 II) . This orbital variability could be under- 
stood in the framework of inverse Compton scattering of photon 
fields in a striped puls ar wind model, which predicts pulsed and 
variable HE emission dPetri & Dubusl 1201 II) . In any case, from 
an observational point of view it is not clear yet if LS I +61 303 
contains an accreting black hole, an accreting neutron star or a 
non-accreting neutron star. 

In addition to LS I +61 303, three binary systems that contain 
a massive star and a compact object and display extended radio 
emission have been clearly detected at VHE: PSR B 1259-63, 
LS 5039 and HESS J0632+057. Thes e four systems have a simi- 
lar sp ectral energy distribution (SED; Dubus 2006: iHinton et al.l 
2009), peaking at MeV-GeV energies, and thus are considered 
gamma-ray binaries. Whereas the nature of the compact object in 
LS I +61 303, LS 5039 and HESS J0632+057 remains unknown, 
PSR B1259-63/ LS 2883 is the only g amma-ray binary with a 
confirmed pulsar (Johnst on et al.ll 19921). T his binary system con- 
tains an 09.5 Ve (iNegueruela et al.ll201 ll) donor and a 47.74ms 
radio pulsar o rbiting it every 3.4 years in a very eccentric orbit 
with e = 0.87 dJohnston et alJll992t Ijohnston et al.l 11994 . The 
radio light curve of the unpulsed emission of PSR B 1259-63 is 
well explained by the adiabatic expansion of a synchrotron bub- 
ble, which is similar to the behaviour found in LS I +61 303. 
Pulsed radio emission is detected with a spectral index of about 
-0.6, although the radio pulses vanish for some weeks dur- 
ing t he periastron passage dConnors et al.l 120021: 1 Johnston et al.l 
2005) probably due to free-free absorption by the stellar wind 
and interaction with the Oe-star disk. Note that in the case of 
LS I +61 303 the apastron separation is smaller than the perias- 
tron separation for PSR B 1259-63, therefore detecting pulses 
from LS I +61 303 is unlikely even if the pulsar exists and 
is beamed at the Earth. The system was detected by HESS in 
the 2 004 and 2007 periastron passages dAharonian et al.l 12005, 
120091) . Extended and variable radio emission at milliarc second 
scales, outside the binary system, was o bserved after the 200 7 
periastron passage of PSR B 1259-63 dMoldon etal.1 l2011al) . 
This is the first observational evidence that a non-accreting pul- 
sar orbiting a massive star can produce a milliarc second ra- 
dio structure similar to those of LS I +61 303, LS 5039 and 
HESS J0632+057. 

As in the case of LS I +61 303, the nature of the compact ob- 
ject in LS 5039 is unknown: it is either a black hole or a neutron 
star with a mass between 1.5 and IOMq in a slightly eccen- 
tric (e = 0.35) 3.9-day orbit around an 06. 5 V((f)) star located at 
2.5 kpc (ICasares et alJl200Sbt IXragona et al.ll2009l) . Similarly to 
LS I +61 303 and PSR B 1259-63, LS 5039 also shows X-rays, 



HE and VHE emission modulated by the o rbit (Takah ashi et al.l 
120091: lAbdo et al.l l2009bl: lAharonian etaD 120061: respectively). 
The discovery of a bipolar extended milliarc second radio emis- 
sion morphology wit h VLBA observatio ns prompted a micro- 
quasar interpretation dParedes et al.ll200Ci|) . but additional VLBA 
observations revealed a changing behaviour of the morphol- 
ogy that cannot e a sily be explaine d by a microquasar scenario 
dRibo et al.l 2008). Rea et al.l (120111) performed recently a deep 
search for pulsations from LS 5039 with Chandra, finding no 
periodic signals in a frequency range of 0.005-175 Hz. 

HESS J0632+057 h as recently j oined t he short list 
of gamma-ray binaries. Mold oriet al.l (1201 lbl) reported the 
slightly extended and variable radio milliarc second structure 
of HESS J0632+057 at 1.6 GHz with the e-EVN. It is po- 
sitionally coincident with a BOpe star and the X-ray source 
XMMU J063259.3+054801 dHinton et al.ll2009l) whic h is also 
variable ( Acciari et al. 2009b). Bonsiorno et al. (2011) have re- 
cently confirmed its binary nature by the discovery of a period- 
icity of 320 + 5 days in a 0.3-10 keV light curve obtained with 
Swift. This faint and point-like VHE s ource was first detecte d 
in the HESS Galactic Plane Survey (lAharonian et al.l 120071) . 
HESS J0632+057 is clearly variable also at VHE: 2006-2007 
and 2008-2009 observations by VERITAS imposed flux upper 
limits well be low the values dete cted by the HESS Galactic 
Plane Survey (lAcciari et al.ll2009bl) . and afterwar ds higher TeV 
gamma-r ay emission wa s detected by VERITAS ( Ong 201 1) and 
MAGIC (lMariottill201ll) . No pulses have been detected at any 
wavelength from HESS J0632+057 dRea & Torresll201ll) . 

Overall, there are many similarities between PSR B 1259-63 
and the other three known gamma-ray binaries (LS 1 +61 303, 
LS 5039 and HESS J0632+057), but no pulsations have been 
detected up to now in the latter. In this paper we report our search 
for radio pulsations from LS I +61 303. In Sect. [2] we present 
the situation of the searches for pulsations from LS I +61 303 
before our work. In Sect. [3] we describe how we selected the 
orbital phase of the observation and the frequency in order to 
minimise the effect of absorption and pulse broadening in the 
observations described in Sect. |4] We describe the data analysis 
in Sect. [3 The results are presented in Sect.0 where we describe 
the implications for future observations. In Sect. [7] we discuss 
these results and present our conclusions. 

2. Previous searches of pulsations 

To our knowledge, only two limited searches for evidence of 
a pulsar in LS I +61 303 were reporte d in refereed journals be- 
fore our observations. ICoe eta! 1 (119821) recorded 10. 7-GHz radio 
data with a time resolution of ~ 15 s, but no evidence of short- 
term periodicities was found in the range 28-4200 s. The only 
ob servational hint of pu lsed emission has been the discovery 
by iPeracaula et al.l (1 19971) of a ~ 1 .4 h periodicity in a series of 
~ 4 mJy microfl ares active during 8 h in a 5 GHz VLA radio ob- 
servation. McS wain et al.l (1201 ll) recently reported on GBT ob- 
servations placing deep upper limits between 4.1 and 14.5 //Jy at 
C, S and X bands. 

LS I +61 303 does not appear as a radio pulsar in any 
of the public catalogues built after performing blind pulsar 
searches. LS I +61 303 has a relatively high dispersion mea- 
sure (DM): the NE2001 electron density model (ICordes & Laziol 
2002) gives DM = 60 + lOpc cm -3 in the direction and distance 
to LS I +61 303 (d = 2.0 + 0.2 kpc). As shown in Appendix[A] 
the pulse broadening is also considerable at low frequencies. 

The Green Bank and Jodrell Bank observatories have per- 
formed pulsar surveys pointing to the LS I +61 303 region. 
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The G reen Bank short-period survey conducted by lStokes et al.l 
(1985) used a sampling time of 2 ms. These authors computed 
a minimum detectable flux density of ~ lOOmJy for a period of 
-4 ms, or ~ 2mJy for P> 100 ms. More recently, the GBT350 
Survey of the Northern Galactic Plane for Radio Pulsars and 
Transients, using a fast sampling time (f sam p = 81.9 2/is), attained 
a ~2 mJy sensitivity for normal, slow pulsars (He ssels et al.l 
2008; this sensitivity should be extendible to 10- 100ms pulsars 
with moderate DMs, J. W. T. Hessels, private comm unication). 

In the case of Jodrell Bank. lClifton & Lynd (fl986) claimed a 
sensitivity of ~ 1 mJy (considering a 6<x detection threshold) for 
a pulsar with a duty cycle of 4% for the Jodrell Bank B survey. 

The surveys of Molonglo, Parkes and Swinburne were made 
from the southern hemisphere, at latitudes where LS I +61 303 
is not visible. LS I +61 303 is outside the range of available 
declinations for Arecibo. 

At other w avelenghts, deep Chandra observations performed 
by lRea etaL l d20ioh did not exhibit X-ray pulses. Despite show- 
ing a HE spectrum similar to those observed in gamma-ray pul- 
sars, Fermi/LAI has not reported pulsed gamma-ray emission 
either. However, the lack of good orbital ephemeris and an a pri- 
ori knowledge on the putative pulsar period hinders the search 
in the currently available Fermi/LAI data. 

3. Minimising absorption 

LS I +61 303 is an eccentric binary system in which the prox- 
imity of the Be star to the compact object should cause a sig- 
nificant free-free absorption of the putative pulsar radio emis- 
sion. This is strongly modulated along the orbit through the pro- 
nounced eccentricity of the system. In our model we consider 
the effect of free-free absorption with the stellar wind of the 
massive star, and we ignore the interaction with the circum- 
stellar disk because it is truncated at 34-37 R , much lower 
than the distance at the considered orbital phase. We consider 
a smooth and continuous stellar wind as a first-order approxi- 
mation although a realistic wind is affected by turbulence and 
may be clumpy. We do not discuss other possible influences of 
the Tsytovitch- Razin effect and the synchrotron self- absorption 
(see for instancelHornbv & Williamsll966tlWhite & Chenl 19951: 
iDoughertv et al .1120031). Other eff ects may be important, as dis- 
cussed in lThompson et al.l (Il994l) . such as refraction of the radio 
beam, pulse smearing, and induced Compton scattering. These 
effects are unknown due to the lack of information about the 
putative pulsar properties and the conditions of the plasma sur- 
rounding it. 

Because we aimed to observe at the orbital phase with a min- 
imum absorption, we computed the optical depth towards the 
compact object at different orbital phases and f or different f re- 
quencies using the Rosseland mean opacity (see Dubus 2006): 



t^0.3 x 
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x 



Jl 



r~ 4 dl . 
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The assumed paramete rs for the Be star are a mass-loss rate of 
M w = 10" 8 M Q yr" 1 (Howarth & Prinia 1989), a terminal wind 
speed o f Vw.oo = 1750 km s" 1 |Hutchings 197JJ) and a /?-law with 
£ = 0.8 (lPulsetal.1 119961) . and a temperature of T = 28 000K 
(ICasares et al.ll2005air Distances are in astronomical units. We 
integrated along the line-of-sight path from the position of the 
compact object (/q) up to the observer at qq . We ha ve consid- 
ered the orbital parameters in Casar es et al ] <2005a). Since we 
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Fig. 1. Expected LS I +61 303 pulsed flux density along the orbit 
at different frequencies for a = -1.0 assuming that the intrinsic 
pulsed flux density at 1.3 GHz is So and considering free-free 
absorption effects. Periastron and apastron are indicated by the 
vertical black lines at phases 0.23 and 0.73, respectively. The 
highest observable flux density is expected to occur at 0.6 GHz 
at orbital phase ~ 0.6-0.7, three days before apastron. The ob- 
servation phase range of this work is indicated by the vertical 
green stripe centred at 0.54. 



are assuming that the compact object is a neutron star, we used 
1.4 Mq for its mass. Hence, the mass function implies an incli- 
nation angle i = 65°, which is slightly out but not very far of 
the rough lim i t of 60 ° imposed by the lack of clear shell lines in 
ICasares et al.l (l2005al) . Note that these assumptions in the stellar 
wind and the choice of the orbital parameters plays an important 
role in the resulting absorption. 

With these opacities we have estimated the absorbed pulsed 
flux densities as a function of the unabsorbed pulsed flux density 
S at 1.3 GHz, 



Soe" T . 



(2) 



at different frequencies and for several values of the spectral in- 
dex (a, defined as S oc y a ): a = -0.5 (a very flat spectrum), -1.0 
(an intermediate value) and -2.0 (a steep spectrum). In Fig.[T]we 
show the expected pulsed flux densities for a = - 1 .0 for different 
frequencies. The modelled flux density variability of the pulsed 
emission at a given frequency is entirely due to the change in 
the the optical depth along the orbit. We conclude that the best 
orbital phase range to try to detect pulsations in LS I +61 303 
is ~ 0.6-0.7, and the highest absorbed flux densities assuming 
a = -1.0 would be received at ~ 0.5-2 GHz. In Fig. [2] we plot the 
opacity for typical values of the stellar parameters. We can see 
that, under the conditions considered here, the opacity is around 
or below 1 at the orbital phase range mentioned above. 



4. Observations 

The Giant Metrewave Radio Telescope (GMRT) is an array of 
30 fully steerable 45 m diameter parabolic dish antennas. It is 
located at 19°6' N, 74°3' E, at a height of 650 m above the sea 
level, around the village of Khodad, 80 km north of Pune (India). 
Fourteen antennas are distributed in a 1000 m central core and 
the remaining are arranged along three ~ 14km-long arms. 
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Fig. 2. Optical depth in logarithmic scale (colour bar) for differen t combinations of ste llar parameters at orbital phase 0.54. The 
dependence with the system inclination, i, in the wide range 15-90° (Casaresetalj2005a) is very low at this orbital ph ase, although it 
is stro ng around periastron. The terminal wind velocity is plotted between 600 km s" 1 (wind veloci t y at the equator in McSwai n et al 



201 1 )) and 2 000 km s" 1 . The expected stellar temperature ranges from 21 500 K lMcSwain et al.l (1201 lb to 30 000 K. iCasares et al 
2005a) estimated the stellar mass to be between 10 and 15 M©. The black circles mark the values adopted in this paper. 



Phased array observations were performed around midnight 
between 2009 July 2 and 3 (UTC), centred at phase 0.54. 
Three observations of LS I +61 303 were carried out, lasting 
171 minutes (21:31:45-23:22:16 UT), 123 minutes (00:52:04- 
03:01:02 UT) and 34 minutes (03:09:35-03:43:06 UT), i.e., a 
total of 5.47 hours. The array was phased on 3C48 twice: the first 
time before any pulsar observation, and the second time between 
the first and the second pointing on LS I +61 303. We decided 
to observe at 1280 MHz to minimise the effect of absorption 
(see Sect. [3]) and avoid the pulse broadening (see Appendix |A| 
and radio frequency interferences (RFI), which are more of- 
ten present at lower frequencies. Because we aimed for a good 
phase coherence between antennas, we did not use the farthest 
three antennas in each arm. Data from two polarisations were 
taken in a bandwidth of 32 MHz centred at 1280 MHz divided 
into two subbands, the upper side band (USB, from 1280 to 
1296MHz) and the lower side band (LSB, 1264 to 1280MHz), 
each split into 256 channels. Three back-ends were used si- 
multaneously: the multi-bit phased array back-end (PA), which 
is specifically conceived for pulsar observations, at USB, the 
polarimeter (PMT) at USB and the PMT at LSB. This redun- 
dancy is convenient because each back-end introduces different 
noise features and this allows cross-checking. The PA sampling 
time was chosen to be t samp = 128 jus, and f samp = 256 /is for both 
PMTs. 

With the aim of checking the performance behaviour and 
data quality, we observed four known pulsars: PSR B0329+54, 
PSR B1937+21, PSR B2022+50, and PSR B2035+36. This 
group of pulsars covers a wide range of the parameters space, 
e.g. periods ranging from 1.56 to 715 ms and fluxes from 0.8 to 
203 mJy. In Table [H obtained from the ATNF Pulsar Catalogued 
(Ma nchester et all2005l) . we show the period CP), the dispersion 
measure (DM), the pulse width at 50% (W 50 ) and 10% (Ww) of 
the flux density peak, the mean flux density at 1400 MHz (S 
and the spectral index. 

We used the four known pulsars to compare the expected 
and the ob t ained S /N in a. manner similar to that described in 
iJoshi et aD (|2009). In Table[2]we show the expected and obtained 
S IN in our observations, as well as the ratios between both val- 



Table 1. Parameters of the four known pulsars from ATNF when 
available. P is the period, DM is the dispersion measure, W50 
and W\q are the width of pulse at 50% and 10% of the peak, 
respectively, is the mean flux density at 1400 MHz and a 
is the spectral index defined as S oc y a . 



Pulsar Name 


P 


DM 


w 50 


w l0 


o mean 
° 1400 


a 




(ms) 


(pc cm -3 ) 


(ms) 


(ms) 


(mJy) 




B0329+54 


714.5 


26.833 


6.6 


31.4 


203 


-1.6 


B1937+21 


1.558 


71.0398 


0.063 


0.19 


10 




B2022+50 


372.6 


33.021 


4.7 


14.9 


2.2 


-0.8 


B2035+36 


618.7 


93.56 


10.7 




0.8 


-1.6 



Table 2. Expected (computed with Eq. ©) over obtained signal- 
to-noise ratios (S/N) for the four known pulsars. 



pulsar name 


(S 1 N) expected 


(S /Attained 


/AOexpected 
(S /^obtained 


PAUSE 








B0329+54 


9241.2 


179.1 


51.60 


B1937+21 


169.81 


31.3 


5.43 


B2022+50 


64.737 


28.9 


2.24 


B2035+36 


32.219 


18.3 


1.76 


PMT USB 








B0329+54 


9618.6 


191.1 


50.33 


B1937+21 


158.33 


49.8 


3.48 


B2022+50 


64.737 


28.8 


2.25 


B2035+36 


32.884 


18.5 


1.78 


PMT LSB 








B0329+54 


8606.1 


167.5 


51.38 


B1937+21 


141.27 


64.3 


2.40 


B2022+50 


59.499 


11.8 


5.04 


B2035+36 


30.181 


12.3 


2.45 



ues, for each of the four known pulsars. The expected S /N were 
computed using 



http ://w w w. atnf.csiro.au/research/pulsar/psrcat 



S/N = 



(3) 



GN a -yjNptAv 
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w 
p ' 



where S mQan is the mean flux density in Jy, (5 is a dimension- 
less constant that accounts for the losses in the digitisation of 
the received analogue signal, r sys = r rece i V er + ^sk y + Aground is 
the total system temperature in Kelvin, G is the single antenna 
gain in K Jy -1 , N a is the number of antennas, N p is the number 
of orthogonal polarisations, t is the integration time in s, Ay is 
the bandwidth in Hz, and D is the duty cycle defined as D - 
where W is the width of the pulses and P their periodicity. 

As explained above, we used N p = 2 and Ay = 16 MHz. We 
considered p = 1/0.8, r sys = 73 K and G = 0.22 KJy" 1 following 
the GMRT System Parameters and Current Statufl document. 
The mean flux density at 1280 MHz was computed using 

and a (S|^ = S^y x (^f) from TableHJwhen a was avail- 
able and assuming a = -1.8 for PSR B 1937+21. We used 20 an- 
tennas at USB and 18 at LSB during the observations before 
3:02 UT (PSR B1937+21, PSR B2022+50, PSR B2035+36 and 
the first two observations on LS I +61 303) and one antenna less 
for the remaining observations (PSR B0329+54 and the last ob- 
servation on LS I +61 303). The duty cycle was computed as 
D = w hen W m was available. For PSR B2035+36 we used 
D = ^ instead. 

The ratio of expected S /N versus obtained S /N should be 
greater than unity because Eq. ^ does not model all error 
sources. This ratio can be much higher than 1 if the data are af- 
fected by corruption due to RFI and deviation from phasing over 
time. In order to take this into account we scaled the sensitivity 
of LS 1+61 303 observations using a parameter that averages the 
expected over obtained S /N ratios, which we will call the cor- 
rection factor K. Hence, the minimum detectable mean pulsed 
flux density is given by 



omean 
°min 



Kx(S/N) n 



sys 



GN,^N^ 



D 



l-D 



(4) 



It can be seen in Table [2] that the LSB data have a lower 
S/N (both expected and obtained) than the USB data. This is 
mainly because two antennas less were used in LSB than in 
USB. Added to this, the PA data have a shorter sampling time 
than PMT (128//S versus 256 yus), we conclude that the PA data 
are the most suitable to detect pulses from LS I +61 303. 

Tabled shows that PSR B0329+54 has a much higher value 
of the expected over obtained ratios compared to the other pul- 
sars. Since PSR B0329+54 is extremely bright, this result is 
probably reflecting a dynamic range problem. We do not expect 
LS I +61 303 to emit extremely bright pulses, therefore we did 
not consider PSR B0329+54 when computing the correction fac- 
tor K. 

For pulsar PSR B 2035 +36 the expected over obtained ra- 
tio is the lowest in each back-end. W\$ is unknown for this 
pulsar; this can lead to a wrong estimation of S /N. Moreover, 
PSR B2035+36 is very faint, near detection threshold, which 
can yield a wrong value of the obtained S IN. We decided not to 
include PSR B2035+36 in the estimation of the correction factor 
K. 

For these reasons we computed K by averaging only over 
PSR B 1937+21 and PSR B2022+50 (see values in Table0). The 
result is K * 3.8 for the PA data and K ^ 2.9 for the PMT USB 
data. 




Fig. 3. Flux density upper limits in mJy in the direction of 
LS I +61 303 obtained with Eq. © for increasing values of the 
duty cycle and for different values of the correction factor K for 
the PA data. 



5. Analysis 

The data of LS 1+61 303 were analysed using the publicly avail- 
able pulsar analysis package SIGPROC0 adapted for the GMRT 
data format. The data were first examined for radio frequency 
interference (RFI) by estimating the distribution of powers over 
the whole time series. Any outliers to the expected Gaussian 
distribution were clipped. The resulting 256-channel data for 
a sideband were dedispersed to 256 trial dispersion measures 
(DM-integrated electron column density in the line of sight) 
ranging from to lOOOpccm -3 . The Fourier transform of the 
dedispersed time series for each trial DM was then computed 
and up to 32 harmonics in the resulting spectrum were summed. 
Any periodicity above eight times the signal-to-noise ratio in the 
summed spectrum was recorded. After ignoring the known inter- 
ference periodicities, the data were folded at candidate periodic- 
ities to produce eight 2-MHz subbands and eight subintegrations 
over the entire observations as well as the total average profile, 
which were plotted as a composite plot. These plots were man- 
ually examined for all candidate periodicities, which were all 
rejected because their DM profiles suggested a terrestrial origin 
and they were not visible in all subintegrations (i.e., they were 
not present during all observations). Hence, no pulsed emission 
from LS I +61 303 was detected, which places a strict upper 
limit on the absorbed pulsed emission from this binary system. 



6. Results and implications 

A pulsed signal with expected DM for LS I +61 303 was not 
detected. Considering K = 3.8, t = 3 hours, a (S /N) min = 8, and a 
canonical duty cycle of 10% in Eq. © we obtained an upper 
limit of 0.38 mJy for the radio pulse mean flux density from a 
putative pulsar with P>2ms in the direction of LS I +61 303 
for the PA USB data from GMRT. In Fig. [3] we show the min- 
imum mean flux density detectable in our observations for dif- 
ferent values of the duty cycle and the correction factor K. For 
the PMT USB data the correction factor is K ^ 2.9 and the cor- 
responding upper limit is 0.27 mJy for pulses with P >5 ms. 

We show in Fig. 4 our obtained upper limit of 0.38 mJy 
at 1.28 GHz and the corresponding upper limits as a function 



http ://w w w.ncra. tifr.res .in/ gtac/GMRT- specs .pdf 



http ://sigproc . sourcef orge . net 
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of frequency for different spectral indexes including absorp- 
tion (solid lines) and without absorption (dashe d lines). In this 
lot w e also show the upper limits obtained by IMcSwain et al.l 
201 lh . even though the observations were not performed at 
the same orbital phase and the sensitivities were computed us- 
ing different methods. We plotted potential ten-hour sensitiv- 
ities for five observatories for future observations planning: 
the Low Frequency Array (LOFAR) considering the HBA; 
the Green Banck Telescope (GBT) using the GBT Ultimate 
Pulsar Processing Instrument (GUPPI) and, for comparison, 
the GBT350 Survey sensitivity; the Effelsberg 1 00-m Radio 
Teles cope considering the surveys presented in iKlein et a"D 
(2004]); the Expanded Very Large Array (EVLA); and the 
Atacama Large Millimeter/submillimeter Array (ALMA). For 
EVLA and ALMA, we considered a putative pulsar back-end 
with a full bandwidth and a reduced bandwidth (Ay = 800 MHz 
as for GUPPI), N p = 2 and ft = 1. We have considered the most 
favourable conditions, and hence these estimates must be con- 
sidered lower limits to the minimum mean flux density of a de- 
tectable pulsar in a ten-hour observation in the LS I +61 303 
direction. The duty cycle is assumed to be D = 10% and in the 
case of the GBT it has been corrected of pulse broadening as 
explained in Appendix lAl 

Looking at Fig. 4 it becomes evident that LOFAR operates 
in a zone where pulses from LS I +61 303 would be totally 
absorbed, and pulse broadening is very important if the puta- 
tive pulsar is fast. In the other extreme of the radio spectrum, a 
ten-hour observation with ALMA would be interesting only in 
the case of a very flat intrinsic spectrum. It would be better to 
observe at the favoured frequency zone from 0.5 to 5 GHz, es- 
pecially in case of a steep spectrum. The Effelsberg telescope 
operates in this frequency range and would provide good sen- 
sitivities. According to Fig. 4, the EVLA would be capable to 
detect the lowest flux densities, but these sensitivities are com- 
puted using many assumptions, considering future planned ca- 
pabilities and an unplanned pulsar back-end with assumed spec- 
ifications. On the other hand, GBT enjoys the good performance 
of GUPPI. Furthermore, the minimum of the GBT curve is at 
lower frequencies, precisely at the few GHz favoured zone. In 
this context, we note that a high-frequency search for pulsations 
from LS I +61 303 has been performed with GBT around orbital 
phase 0.9, where absorption effects are not important at these 
frequencies according to our model, placing de ep upper limits 
betwe en 4.1 and 14.5 //Jy at bands S, C and X (IMcSwain et al.l 
I20TH) . 



7. Discussion and conclusions 

Detection of pulses is the only definitive way to confirm the pul- 
sar nature of the compact object in LS I +61 303. We have per- 
formed phased array observations with the GMRT at 1280 MHz 
in the LS I +61 303 direction centred at phase 0.54, from 
21:31:45 2009 July 2 UT to 03:43:06 2009 July 3 UT. No pulses 
have been found in the data set, with a minimum detectable flux 
density of around 0.38 mJy for periods greater than 2 ms and 
a duty cycle of 10%. If LS I +61 303 had the same duty cy- 
cle as PSR B 1259-63 (63%), the upper limit would be 1.3 mJy 
(see Fig. [3]), which is belo w the ~ 5 mJy mean flux density of 
this pulsar at this frequency (I Johnston et al. 112005b . These results 
give a guide for future observations planning, which would be 
best at ~ 0.5-5 GHz because the upper limit at 1.28 GHz implies 
deep upper limits at high frequencies in case of a steep spec- 
trum, as shown in Fig. 4. Following the discussion from Sect.[3j 
the orbital phases at which the flux is higher are ~ 0.6-0.7, and 



at these phases the higher flux densities are received at -0.5- 

2 GHz assuming a spectral index a = -1.0. In addition to this, 
below 0.5 GHz the pulse broadening due to the ISM scattering 
is very important for fast pulsars (see Appendix lAl for more de- 
tails). Therefore, long (~ ten hours) observations at ~ 0.5-5 GHz 
with Effelsberg, EVLA or more likely GBT directly pointing to 
LS I +61 303 around phases 0.6-0.7 are recommended to signif- 
icantly improv e these upper limits or have a good chance to de- 
tect the pulses. McS wain et al.l (l201ll) performed a pulsar search 
in LS I +61 303 using GBT at a similar frequency range (bands 
C, S and X) and a different orbital phase (0.9), placing deep up- 
per limits (4.1-14.5/Jy). 

In this argumentation we are not considering the possibil- 
ity of the pulsar beam not pointing close to our line of sight. 
This issue is less problematic at X-rays, because the beam is 
usually wider than the radio one. However, no X-ray pulsations 
have been detected so far from LS I +61 303. A deep search 
performed with Chandra observations provided 3<x upper lim- 
its on the pulsed fraction of the total X-ray emissio n ranging 
betwe en 7 and 15% in the 0.3-10 keV energy band dRea et al.l 
120101) . It should be mentioned that the TeV binary system host- 
ing PSR B 1259- 63 does not show X-ray pu lsations during its pe- 
riastron passage (iChernyakova et al. 2006), with an upper limit 
of 15% on the pulsed fraction (IRea et al. 2010), yet it is known 
to be a pulsar beamed at Earth. 

An alternative way to experimentally prove the presence of 
a pulsar in LS I +61 303 is finding HE pulsations. As mentioned 
in Sect. [T] the HE emission detected by FermifLAI could be 
explained by inverse Compton scattering of photon fields in a 
striped pulsar wind model, which predicts pulsed and variable 
HE emission (Petri & Dubusll201lh . 

A better orbital solution would improve the opacity estima- 
tion and the knowledge of the acceleration suffered by the puta- 
tive pulsar and therefore the chances of resolving the pulses in 
observations along a large part of the orbit. In any case, because 
of the opacity of the binary system and the possibility of the pu- 
tative pulsar not beaming towards Earth, it might be impossible 
to detect pulses even with infinite sensitivity. 
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Appendix A: Pulse width and broadening 

The ISM electron density variations cause a multipath scatter- 
ing reflected in a temporal broadening of the pulses at high DMs 
and low frequencies ( see for instance [C ordes & Rickett 1998; 
ICordes & Lazioll2002l:lLambert & Ri ckett 1999). In our case we 
have corrected the duty cycle in Eqs. © and © using the fol- 
lowing method. We have computed the flux density as a function 
of time, S(t), convolving an intrinsic pulse profile, Si nt (t), and a 
broadening function b(t). Only with this purpose we considered 
an intrinsic Gaussian pulse profile with standard deviation <x, 



p^ mean 
Sint(0 = e 



V2 



(A.1) 



with a period of P = 10 ms, an intrinsic duty cycle of Ant = 10% 
and hence an intrinsic pulse width of W[ nt = 1 ms. In order to link 
o~ and Wint in this context, we considered W - Wio (the pulse 
width at 10% of the peak). Therefore, cr = Wi 0Mt / V81nl0. 

The temporal pulse broadening is parametrised with the 
timescale Td, defined as the first moment of the broadening func- 
tion, roughly a one-sided exponential, determined by deconvolv- 
ing the intrinsic and measured pulse shapes ( see for instance 
ICordes & Laziol 120021: ICordes & Ricked 1 1998b . We effectively 
considered a one-sided exponential for the pulse broadening 
function b(t): 



b(t) = 



^-e T d 




t >0 
t<0. 



(A.2) 



Hence, Td is the width of the pulse to 1/e of the maximum 
pulse amplitude if the intrinsic pulse shape is a delta function. 
We have obtained the values of Td at several frequencies u sing 
the NE2001 electron density model (ICordes & Laziol 120021) in 
the direction of LS I +61 303 and considering a distance to the 
source of d = 2.0 ± 0.2 kpc. 

Finally, the measured pulse profile S (t) is the convolution of 
Smt(t)andb(t): 

S(t) = S int (t)®b(t) = 



2(TT d 



x e 



x erfc - 



(T 



(A.3) 



V2o- V2T d , 

In Fig. IA.1I we show the resulting pulse widths with 
P = 10 ms and 50 ms, and Ant = 10% (top panel) and 50% (bot- 
tom panel). Vertical lines indicate the frequencies v c at which 
D = 75% and thus the pulse is hardly detectable. For an intrin- 
sic duty cycle of 10%, above ~ 0.5 GHz the effect is negligible. 
We conclude that the pulse broadening does not affect our ob- 
servations at 1.28 GHz. In Fig. 4 this correction is used consid- 
ering P= 10 ms to the points below 1 GHz. According to this 
result, at 120 MHz (the High Band of LOFAR) this presumed 
pulsar would not be detectable because the pulse width would 
be greater than the period. For this reason the LOFAR point 
is not corrected for pulse broadening. If the compact object in 
LS I +61 303 is a fast (P~ 10 ms) pulsar, observations above 
0.5 GHz are recommended to avoid the pulse broadening by the 
ISM. 
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Fig. A.l. Pulse widths broadened due to ISM scattering in the 
frequency range at which the role of pulse broadening turns from 
crucial to negligible. The resulting width W\q is plotted consid- 
ering J P=10ms (blue dashed lines) and ,P=50ms (green dash- 
dotted lines) for intrinsic duty cycles Ant =10% (top panel) and 
Ant =50% (bottom panel). The pulse broadening timescale Td is 
plotted with a solid indigo line. Horizontal lines mark 75% of 
period, and vertical lines indicate the v c where the pulse widths 
cross this limit. 
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